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ABSTRACT 
The thermal structure of the ocean, especially the uppermost mixed 
layer, greatly affects sonar ranges. In this paper, similarity theory 
is applied to the problem of forecasting the depth of the mixed layer 
during the warm season, assuming the controlling processes are secular, 


non-advective, and non-divergent. The resulting forecast method consists 
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mainly of two equations. Parameters used are wind, coriolis effect, the 
coefficient of thermal expansion and a measure of the excess heat within 
the mixed layer. The constants in the equations were determined using 
data from OWS Papa (50N, 145W). The forecast method treats both seasonal 
and transitional thermoclines. The method was tested with data from OWS 
Papa and OWS November (30N, 140W). The tests apparently indicate wide 
applicability of this forecast method and thus tend to corroborate the 
proposal by Kitaigorodsky that the mixed-layer depth is a function of a 


universal coefficient. 
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l. introduction. 

The detection of enemy submarines by sonar is one of the major un- 
solved problems of this decade. Despite recent technological advances 
in many fields of science, the nuclear submarine is still practically 
invulnerable. Even though sonar technology, in particular, has advanced 
rapidly, we are still unable to position our ASW forces or their sonars 
for optimum performance. This is Largely because sonar ranges are either 
enhanced or reduced by refraction of sound energy in the surface layers 
and we are not able to forecast accurately the thermal structure which 
greatly determines refraction. 

Many authors have devised systems for forecasting the thermal struc- 
ture of the upper layer of the ocean based on either dynamical analysis 
or empirical relationships. The equations resulting from dynamical 
analysis are either too complicated, if all the physical processes are 
considered; or are impractical, if many simplifying assumptions are made. 
Empirical relations have been determined for certain locations and ger 
limited time periods but do not appear to be valid universally. 

Another possible method for developing a forecasting system uses 
similarity theory. Kitaigorodsky [1] investigated the application of 
similarity theory to the ocean thermal-structure-forecasting problem. 
Although his results have certain drawbacks, as noted in appendix IV, 
still he has shown the applicability of the method. Consequently, the 
method used by Kitaigorodsky is also applied by this author, with some 
modifications of parameters, in an effort to develop a more practical 
result. The form of P, a dimensionless coefficient inherent in the 


application of similarity theory, is determined by use of data from 








Ocean Weather Ship Papa (50N 145W). 

Figure 1 depicts a typical ocean thermal structure in the warm sea- 
son at OWS Papa. Characteristically, there is a guasi-isothermal layer 
that extends from the surface to the upper boundary of a negative tem- 
perature gradient. 

The depth to which the quasi-isothermal layer extends is usually 
referred to as the mixed-layer depth. In this paper the mixed-layer 
depth (MLD) is defined as that depth at which the temperature of the 
water first becomes ic less than that of the water at the surface. 

A transitional layer of negative temperature gradient between layers 
of relatively small temperature variation is referred to as a thermocline. 
Four main types of thermocline can be classified, primarily according to 
degree of permanence: diurnal, transitional, seasonal, and permanent. 

Diurnal thermocline: it results from a net heat gain during 
the day, a small thermocline (AT < 1C, see fig. 1) being formed 
close to the surface by late afternoon. With a net heat loss at 
night, the thermocline will be destroyed by morning. 

Transitional thermocline: a moderately large thermocline 
(AT > IC) is formed when diurnal heat input exceeds losses. After 
a few days or weeks it joins with the seasonal thermocline as the 
added heat diffuses downward. 

Seasonal thermocline: the transition zone which lies between 
the surface waters warmed during the summer and the colder water 
below. At OWS Papa the seasonal thermocline depth is about 20 
meters in the late summer, and about 50 meters in the spring and 


fall. It is not present ín the winter. 
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the Ocean During the Warm Season at OWS Papa 
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Permanent thermocline; in many localities a deep thermocline 
that is relatively persistent in depth exists throughout the year. 
Below a permanent thermocline the temperature gradient (usually 
negative) is very small all the way to the bottom, 

The entire thermal structure of the ocean, though of interest, is 
not immediately as important to sound propagation with present sonars as 
the depth of the uppermost mixed-layer of the ocean. As far as this 
mixed layer is concerned, there are two basic seasons at OWS Papa. These 
are the warm season, when mixing of the upper layer ís mainly due to the 
wind, and the cool season, when mixing is mainly due to thermohaline con- 
vection, Thus the physical processes during the two seasons are quite 
different and are usually treated separately. 

As in the paper by Kitaigorodsky, similarity theory will be applied 
to develop a method for forecasting the mixed-layer depths associated 
with the transitional and seasonal thermoclines during the warm season 
at OWS Papa. Henceforth, diurnal thermoclines will be ignored and the 
terms "MLD" and "depth of the thermocline" will be synonymous and will 
refer only to mixed-layer depths associated with either transitional or 
seasonal thermoclines. Only those changes in the MLD which are secular, 


non-advective, and non-divergent will be considered. 








2. Description of Similarity Theory. 

There are three basic steps in the application of similarity 
theory. 

1) Determine the physical processes that control the physical 

phenomenon of interest. 

2) Select parameters that accurately represent the controlling 

physical processes. 

3) Apply the f theorem [ 2 | to the chosen parameters. 

The 4r theorem is a method for determining a dimensionally-correct 
relationship for a given set of parameters. A short descríption given 
by Binder [3 ) follows. 

Let A,A2A3....A, be n physical quantities which are involved in 
some physical phenomenon. Examples of these physical quantities 
are velocity, viscosity, and density. Let m be the number of all 
the primary or fundamental units (such as length, mass, and time) 
involved in this group of physical quantities. The physical equa- 


tion, or the functional relation between these quantities, can 
be written as 


S(A, A, Agere Am) =O 
The % theorem states that the foregoing relation can be written as 
0) ۱112 رک اا 2 ددم‎ © 


where each W is an independent dimensionless product of some of 
the A's. 


Thus, if there are n physical variables in a particular problem, and m 
fundamental units, then the physical relationship can be expressed in a 
form ivo an) dimensionless ratios. An elaborate and formal 
proof of the theorem is given by Buckingham I. A simplified step-by- 
step procedure and an example are given ín appendix I. 

As can be seen in appendix I, to equate these (n-m) ratios, the 1T 
theorem introduces a dimensionless coefficient (P) which serves to make 
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the dimensionless ratios numerically equivalent. A practical example 
for the necessity of P is given in the second paragraph of section 5. 
The form of P is normally determined from actual test data. P can be a 
constant or it can be a function of some or all of the parameters used 
to determine the ratios. Buckingham [2] suggests several ways to deter- 
mine the form of P. 

One method is to let P be a function of all of the parameters of 
the process with the exception of the one for which the relationship is 
being developed. For ease of reference, this combination of parameters 
is called N. The Tl theorem is entered with the chosen parameters and 
an expression for the dimensionless ratio N is obtained. Pairs of 
values of P and N, determined using observed data, are plotted together. 
Then the form of P(N) can be determined from this plot. 

Once the dimensionless ratios and the form of P have been deter- 


mined, then the results can be used to forecast the mixed-layer depth. 
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3. Variables and Processes that Affect the Mixed-Layer Depth During 
the Warm Season. 
There is general agreement among many researchers that during the 
warming season the mixed-layer depth varies with: 
1) wind (Geary [4], Tabata [5], Mazeika (6], Munk and Anderson [7Ì); 
2) the heat flow between upper and lower layers (Mazeika [6], 
Rossby and Montgomery [8], Kitaigorodsky Ci): 
3) coriolis effect (Rossby and Montgomery (8], Kitaigorodsky CU); 
4) divergence in the upper layers of the ocean (Mazeika L6], Tabata 
and Giovando [9]); 
5) internal wave action (Mazeika [6], Tully [10], Tabata and Gio- 
vando [9]); 
6) advection (Mazeika [6], Tully [10], Tabata and Giovando [9] ). 
Only the first three of these processes will be considered in developing 
a forecast method. The remainder of the influences contribute to the 


scatter of the results. 





4. Selection of Parameters to Represent Controlling Processes. 

An aid in determining representative parameters is to write the 
relevant dynamical equations describing the processes and to examine 
them term by term. 

The equations of motion in the xy-plane, assuming an unaccelerated 


flow in an unbounded and horizontally-homogeneous ocean, reduce to: 


و 0 
کو مر سس سر ه 6 
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(1) 
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where u and v are the components of velocity, -f. 


>= and = are the com- 


ponents of stress on the horizontal plane, and 3 is the coriolis param- 


eter. At the interface, the stress is given by: 
e — e 


where | is the stress on the water surface, W is the mean wind speed and 

Y is a coefficient that is a function of wind speed and the height at 
which the wind is measured. This indicates that the water motions, in- 
cluding che eddies that carry both heat and momentum, depend upon the 
surface wind speed. Thus a measure of wind speed that is representative 
of the stress effects should be introduced as a parameter in the applica- 
tion of the “]| theorem. However, the depth to which wind effects reach 
depends also on the coriolis parameter as shown theoretically by 
Ekman rite Thus, the coriolis parameter must also be considered in 
application of TT theory. 

The maintenance of a thermal structure involves heat flow. A common 


equation for temperature changes due to heat flow is: 
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where /° is the density of water, a is the specific heat of water, T is 
the temperature of the water, and K is the eddy-conductivity coefficient 
for heat. This indicates that the local rate of change of temperature 

is a function of K as well as of the differences in temperature between 
layers. But K depends upon the field of motion and the stability. Param- 
eters to represent motion have already been chosen in the previous para- 
graph. Thus, to complete the representation of the variables in the heat 
flow process, the only additional parameters that need be introduced are 
those representative of stability. 

Thus it appears that, for the conditions and assumptions given in 
section 3, the Tf theorem should be entered with parameters that are 
representative of: 

1) wind stress; 

2) coriolis effect; 

3) stability. 

WIND STRESS. Here it is necessary to find a value of wind speed 
which is representative of the influence of wind on the MLD over a period 
of time. The wind enters into the eas Dyer -depth problem in the form 
of stress, and the stress coefficient (8) is an increasing function of 
wind speed. Thus equations (1) and (2) show that higher winds should be 
weighted more heavily in estimating the influence of wind stress upon 
the MLD. 

One therefore should avoid a linearly-averaged wind. The author 


has devised an empirical method which qualitatively and objectively 





weights the stronger winds in order to obtain a suitably representative 
wind parameter. The parameter is called the "representative maximum 
wind" and is defined as: 
the average of the five highest winds of the eight usually 

reported during a chosen 24-hour period. The chosen period is 

that one having the highest winds in the interval from 72 hours 

to 12 hours prior to observation time. 
The step-by-step procedure for determining the 'representative maximum 
wind" is given in appendix III. 

CORIOLIS EFFECT. The standard equation for the coriolis parameter 
times 10^ will be used to represent coriolis effects: 


< : 
SING X10‏ لیا 2 = ھا دہ مر 


The constant, 1x 10%, was introduced arbitrarily for convenience. 

STABILITY. The density gradient, and thus the stability, can be 
represented by a combination of the coefficient of thermal expansion 
(3 ) and the temperature difference between layers; the following two 
parameters have been chosen to represent stability: 

1) the coefficient of thermal expansion (/3); 

2) the excess heat (Q) in the mixec layer over that in lower layers. 

Values of /3 used in this paper were interpolated from values given 
by Sverdrup (12) assuming a constant salinity of 32.5 %o, These values 
are tabulated in appendix III. Constant salinity was assumed at OWS 
Papa during the warm season as in Tabata [5]. 

The step-by-step methods used to determine Q are given in appendix 


III. One difficulty with the use of Q as a measure of stability is 


illustrated in figure 2. This shows that thermal structures of obviously 
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Figure 2 


Two Possible Thermal Structures with. 
Same Excess Heat Present in Warm Layer 
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different stability can prevail with the same Q. However, the combina- 
tion of 0, 2, and the MLD give a good approximation of the stability 
of the mixed layer. (The MLD enters the forecast model in the form of 
P(N) since P(N) is initially determined from raw data that includes 
MIDs. ) 

Thus, besides the MLD itself, four parameters have been chosen for 
entry into the T^ theorem: 

1) the representative maximum wind (W); 

2) the coriolis parameter ($); 

3) the coefficient of thermal expansion (/2); 


4) the excess heat in the upper layer (Q). 
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5. Practical Application of Similarity Theory. 

In section £, four parameters were chosen as representative of the 
processes effective in forming and maintaining a mixed layer at OWS Papa 
in the warm season, When the Tf” theorem is applied to these four param- 


eters the results are: 


MLD = P( a ne (3) 
= Q Ll 
N= J (4) 


Both of these equations are derived in the example in appendix I. 
Because the relationship between MLD and the other parameters is 
complex, the coefficient P(N) is not in general a constant and must be 
found as a function of some of the parameters, as in appendix I. m an 
example of the ambiguity which must be eliminated by use of P, consider 
again the parameter 0; as shown in fig. 2, two guite different thermal 
structures may have the same excess heat Q. Thus a value of Q alone 
does not uniquely specify the thermal structure. Rather the structure 
is a result of the interaction of all the controlling parameters. That 
is, a wide range of — mixed-layer depths exists for each value of 
a particular controlling parameter, whereas only one MLD exists for each 
combination of parameters. To specify the form of P(N), which is a 
function of all the prevailing parameters, several pairs of values of 
P and N are determined from experimental data and plotted together. 
Then the form of P(N) is given by this plot. 
If,in the application of similarity theory, parameters truly repre- 


sentative of the controlling processes are chosen, then the plot of P 
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versus N should have little scatter. If the parameters used are not 
truly representative, then large scatter results. 

To test equations (3) and (4) they were used to obtain paired values 
of P and N; these values were based on mean monthly data at OWS Papa 
tabulated by Tabata and Giovando [9]. Values of the parameters, W, 7 و‎ 
LL, and Qo (where the subscript S refers to seasonal thermoclines) were 
determined in accordance with Appendix III. The paired values of P and 
N are plotted in figure 3. The points have very little scatter; and, in 
fact, a straight line gives a very good least-squares fit. However, the 
use of figure 3 is limited because it represents a climatological MLD, 
whereas the MLD for a particular day can vary considerably from the clima- 
tological MLD. Figure 3 is important in that the small scatter of the 
paired values apparently indicates that the proper parameters have been 
used in equations (3) and (4). 

In order to make equations (3) and (4) useful in short range predic- 
tion of the MLD, P(N) was determined using daily data from OWS Papa. The 
data were from June to October during the years 1958 through 1962. In 
order to filter out internal wave effects, data were used only from those 
days for which six or more bathythermographs were available; the MLD 
used is a mean of the six (or more) bathythermographs. This requirement 
limited the available data such that only 22 pairs of P and N are avail- 
able for transitional MLDs and 29 pairs for seasonal MLDs. Although the 
number of paired values appears small, they represent over 200 BTs. 

Pairs of values of P and N were first obtained for MLDs associated 
with transitional thermoclines with values of the parameters W, A, 


SVL, and Qr determined as stated in Appendix III. Note that for a 
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P versus N Based on Monthly Climatological Data 
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transitional thermocline the general parameter Q in equations (3) and 
(4) is represented by the more precise parameter Cr. Figure 4 is a plot 
of pairs of P and N resulting from data listed in table 1. The least- 
Squares best-fit straight line has the equation: 


UN 


Substituting this into equation (3) gives: 


MLD =-.25X 10 dou + 2.1 5 Ww) (5) 


Thus equation (5) can be used to forecast a transitional MLD without 
determining a value of P. 

It is interesting to note the similarity between equation (5) and 
the equation derived by Rossby and Montgomery L8]: 


| 6 
MLD=-K,w EKI zT (6) 





W, = wind speed measured at the sea surface. 
Ó = latitude 

Based on 39 observations, they arrived at a value for K, of 2.38 per 
second. It was concluded that K, was negligible. The MLD data for these 
39 observations were obtained by interpolating between readings of tem- 
perature and salinity which are 10-25 meters apart. The average MLD for 
the 39 observations was 32 meters,and they considered a five-meter aver- 
age error as a plausible assumption. Since the data used by Rossby and 
Montgomery appeared to be of poor quality compared to the data available 


at OWS Papa, it is suggested that equation (5) is more applicable than 


(6) for forecasting the transitional thermocline. Tests in section 6 
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FOR 


DATE 


61758 


61858 


62058 


71758 


81658 


71062 


71262 


80661 


90661 
91461 
62362 
62762 
63062 
70262 
70662 
70762 
71462 
71662 
72062 
72662 
72862 


73062 


TASLE 1 


DATA USED TO DETERMINE VALUES OF P AND N 


W 


19.6 


30.6 


12.3. 


26.0 
24.0 
20.0 
25.0 
25.0 
20.0 
20.0 
772 
1372 
17.0 


22.8 


(KNOTS) (Kg cal/cm?) (METERS) 


Or MLD 
7.21 26.4 
6.39 23.4 
8.06 28,5 
3.78 20.4 
18.29 34.1 
9.00 22 
75 30.0 
13.50 23.9 
210 34.8 
24.00 39.7 
2.98 13.9 
5.86 27.3 
4.01 27.4 
6.32 32.4 
9.63 30.3 
10.11 32.4 
5,38 27.6 
5.87 28 .0 
4.86 24.3 
1.10 14.1 
1.83 19.8 
5.34 23.0 
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TS 
(°c) 


10.3 
10.3 
10.4 
13.2 
1225 
10.7 
10.8 


13,5 


95 
10.0 
10.1 
10.4 


10.3 


13.0 
13.0 


W 


P x 10^ 
2.0 
1.6 
2.4 
1.8 
8.6 
4.1 
5.0 
8.0 
SE 


9.4 


3.3 
3.4 


0.8 


MIXED-LAYER DEPTHS ASSOCIATED WITH TRANSITIONAL THERMOCLINES 


Nx 10^ 
0.96 
0.85 
1.10 
0.80 
2295 
1.3. 
1-7 
3.00 
3.40 
3.33 
0.92 
0.86 
6.10 
1.15 
1.40 
1.47 
1.04 
1.13 
1.12 
0.35 


0.46 


E 


- (m o nn 


ho 
© 


= 


" 
Pi 





indicate this to be true and also show that the terms of equation (5) 
are all of the same order of magnitude. 

Although equation (5) apparently is applicable for MLDs associated 
with the transitional thermocline, tests indicated that (5) was no 
longer applicable when the transitional MLD joins the seasonal MLD. In 
this case, a slightly different approach is necessary, apparently due to 
the increased stability associated with a seasonal type thermocline. To 
take this increased stability into consideration Qs (defined ín A»pendix 
III) is used in equations (3) and (4) instead of Q4. 

The general procedure for finding values of P and N for the sea- 
sonal MLD was the same as for values associated with the transitional 
MLD with the exception that Qs was used instead of ۹۰ Pairs of values 
of P and N for seasonal MLDs, based on data listed in table 2, are 
plotted in figure 5. The least-squares best-fit line for these points 


Was; 


p= —61x10 7+ 3.89 N 


Substituting into equation (3) gives: 
4 `. 
MuD = -6.1x10 (ssa) 3.8306) (7) 


Equation (7) is applicable to MLDs associated with the seasonal thermo- 


cline. 

An objective method has been devised to determine whether a transi- 
tional or seasonal condition exists. When N < 3.5 x 1O the thermo- 
cline is defined as transitional. When N > 3.5 x 1075, a seasonal 


thermocline is said to exist. This is not a rigid rule. Obviously, even 
when N «3.5 x 107^, if the forecast MLD usíng equation (5) is deeper 
than a prevailing seasonal thermocline, then equation (7) should be used 


rather than (5}. 
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TABLE 2 


DATA USED TO DETERMINE VALUES OF P AND N FOR MIXED-LAYER DEPTHS 


DATE 


62762 
91958 
61858 
91161 
82761 
82361 
80661 


72062 
90661 
91461 
73062 
90961 
71462 
71662 
91461 
81261 
71262 
71062 
70262 
80861 
81461 
63062 
70662 
70762 
61758 
61958 
62058 
81658 
91858 


W 


24.0 


50.0 
27.4 


27.2 
21.4 
2120 
19.6 


17.2 
26.6 
2702 
23.0 
27.2 
20.0 
20.0 
30.6 
28.4 
21.0 
20.0 
20.0 
28.0 
32.0 
24.0 
25.0 
25.0 
27.4 
27.4 
27.8 
25.6 
50.0 


Qs 2 MLD 
(KNOTS) (Kg cal/cm”) (METERS) 

28.3 27.3 
44.1 49.8 
22 23.4 
43.2 40.5 
40.6 32:7 
41.4 31.3 
35.7 23.9 
28.5 24.3 
39.7 34.8 
40.6 39027 
31.8 22.8 
40.7 37.6 
31.0 27.6 
275 28.0 
39.5 39.7 
3205 30.3 
29.6 30.0 
2155 27.2 
26.2 32.4 
29.6 26.8 
37.4 34.7 
29.1 27.4 
30.5 30.3 
27.6 32.4 
29.5 26.4 
29.8 2525 
29.3 28.5 
39.8 34.1 
43.7 49.4 


ASSOCTATED WITH SEASONAL THERMOCLINES 
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TS 
(°C) 
9.5 


11.6 
10.3 


1.3.3 
14.5 
13.8 
13.5 
11.6 
13.3 
13.4 
12.4 
13.4 
11.0 
11.2 
13.4 
13.2 
10.7 
10.8 
10.1 
13:3 
12.8 
10.0 
10.4 
10.3 
10.4 
10.4 
10.3 
1249 
11.5 


P x 10^ 


19.3 


7.6 
5.6 


21.8 
Zu, 
27.7 
21.0 


20.2 
18.0 
20.1 
123 
LET 
18.0 
16.2 
15.5 
11.3 
16.6 
15.4 
16.9 
59 
11-4 
11.0 
1l. 7 
11.4 
8.2 
8.0 
8.6 
18.6 
7.4 


N x 10^ 


5.3 
5.3% 
49V 
5.3 
5.2 
4.8 
4.5 
4.8 
4.4: 
4.5“ 
4.0 ~ 
3 ٭عۃ‎ 
4.0 ۳ 
3.8 
6.4- 
3.4 
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6. Application and Test of Equations. 

Equations (5) and (7) can be used to forecast MLD over any length 
of time for which the parameters can be accurately predicted. Thirteen 
tests of the equations were made using independent data for forecast 
periods of from one to four days. Five of the tests were based on data 
from OWS Papa (50N, 145W) and eight on data from OWS November (30N, 140W). 
Three of the tests were for transitional thermoclines and ten were for 
seasonal thermoclines. The tests were made using a mean MLD computed 
from four to six BTs from each day as a representative MLD for that day. 
Results are shown in table 3. 

The verification of a forecast is made difficult by the fact that 
the MLD varies randomly. One objective way to measure verification suc- 
cess is to determine if the forecast value lies within the normal range 
of variability of the MLD as measured by the standard deviations of the 
individual values. In table 3 the 13 values listed under "'verified MLD" 
actually represent 63 BTs. The standard deviation with respect to the 
daily mean for these 63 bathythermographs was 3.4/7 meters. The percent- 
age of forecasts that fell within one and two standard deviations of the 
individual MLDs is given in table 4. The root-mean square and the alge- 
braic mean of the forecast differences are also given in table 4, 

The algebraic mean of the forecast was computed to determine if the 
forecasts were biased. Forecasts using the equation proposed by Rcssby 
and Montgomery were biased positively whereas those based on equations 
(5) and (7) were less bíased and negative. This suggests that the nega- 
tive term neglected by Rossby and Montgomery (see section 5) should be 


included as ín equations (5) and (7). 
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TABLE 4 


STATISTICAL RESULTS OF FORECASTS 


% FORECASTS 
WITHIN 1 STANDARD 
DEVIATION 


% FORECASTS 
WITHIN 2 STANDARD 
DEVIATIONS 


ROOT MEAN SQUARE 
OF DIFFERENCES 
(METERS) 


ALGEBRAIC MEAN 
OF DIFFERENCES 
(METERS) 





EI A A مو‎ 


FORECASTS BY 
EQUATIONS 
(5) AND (7) 


=- F ammo A Pb کس‎ 


FORECASTS BY 
PERSISTENCE 


TITO: IESO. EI AA m né ہے کہہے‎ 


69 


92 


24 


46 


6.3 


or -1.06 


tn nenn‏ عص وای 


FORECASTS BY 
ROSSBY AND 
MONTCOMERY 

EQUATION 


u ana An nn 
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To test how the function P(N) at OWS November compares with that 
found at OWS Papa, eight paired values of P and N at OWS November were 
computed. The BTs upon which the verifications were based in table 3 
were also used to obtain these paired values in the same manner as for 
OWS Papa (section 5). Figure 6 is a plot of paired values of P and N 
for the data listed in table 5; the straight line represents P(N) for 
the constants derived from data at OWS Papa. The function P(N) appears 
to be very nearly the same at both OWS Papa and OWS November. 

Two examples of the application of equations (5) and (7) to fore- 
casting follow. 

Example 1. Forecast to be made for 23 August 1957 based on BTs 

taken on 19 August at OWS Papa. 

l. From the BTs taken on 19 August 1957, Gr was determined to be 

1.74 Kg len by methods described in A»pendix III. 

2. Using the August climatological data computed by Tabata and 

Giovando (9], the average net heat transfer downward across the air- 

sea boundary at OWS Papa is approximately +0.1 Kg cal/em? per day. 

Using this information, 0.4 Kg cal/cm* were added to Qr during the 

four-day interval up to the verification day giving: 

Qr = 1.74 + 0.4 = 2.14 Kg cal/cm* 

(Where enough current meteorological data are available, the methods 

for computing the net heat transfer proposed by Laevastu )13[ could 

be applied instead of climatology. ) 

3. Values of /2, NA., and W were determined by the methods described 

in Appendix III. 


4. Equation (4) was entered and N was found to be less than 
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TABLE 5 


DATA USED TO DETERMINE WALUES OF P AND N FOR MIXED LAYER DEPTHS 
ASSOCIATED WITH SEASONAL THERMOCLINES AT OWS NOVEMBER 


DATE MLD W Qs Px 104 N x10 
(METERS) (KTS) ۰ (Kg cal/em“) 
61657 27.9 13.2 32.9 26.8 7.5 
70757 43.0 17.1 31.6 13.0 4.7 
71057 37.5 19.3 22 7 11.5 4.3 
72157 44.5 23.5 120 000 24.9 6.4 
72757 49.0 17.8 28.3 28.5 9.7 
82558 34.7 15.1 39.6 34.2 10.5 
73159 43.5 18.5 33.0 23.1 6.9 
70159 32.0 18.0 18.2 9.9 3.9 
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395 x ose. Conseguently, eguation (5) was used, giving a forecast 
MLD of 16.1 meters. 

5. The seasonal MLD was at 32 meters on 19 August 1957. According 
to the above forecast, the wind was not strong enough to cause the 
transitional thermocline to merge with the seasonal thermocline. 
Thus persistence is used to forecast the depth of this seasonal 
thermocline. 

6. The forecast thermal structure for 23 August 1957 consists of a 
mixed layer down to a transitional thermocline at 16.1 meters with 

a seasonal thermocline below it at 32 meters. The verification 
showed a transitional thermocline at 18.3 meters and a seasonal 
thermocline at 30.5 meters. 

Example 2, For the forecast MLD of 1 August 1959 at OWS Papa, steps 
1, 2 and 3 as done in example 1 resulted in N >3.5 x 09 Conse- 
quently, the procedure for a seasonal thermocline was follo: ed, 

That is, Qs was determined and equation (7) was used to forecast 


the MLD. The forecast MLD was 35 meters and 39 meters verifíed. 
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7, Conclusions. 

The agreement between forecast and cbserved MLDs tends to verify 
the findings of Kitaigorodsky that similarity theory is useful in this 
forecasting problem. Equations (5) and (7) were successful in forecast- 
ing both large and small changes in the MLD, in the positive as well as 
in the negative sense. The concept of a universal function P(N), as pro- 
posed by Kitaigorodsky, is strengthened. 

Equations (5) and (7) appear to be useful in forecasting the mixed 
layer depth to a reasonable degree of accuracy. Some of the deviations 
can be attributed to the effects of divergence, advection and internal 
waves which were not evaluated in this study. 

Future res. arch could well be applied to determining more paired 
values of P and N for other ocean locations. This would not only help 
to fix the constants of equations (5) and (7) better, but would help to 


demonstrate further the universality of P(N), if it exists. 
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APPENDIX I 
BASIC PROCEDURES FOR APPLICATION OF TT THEOREM 

Let n represent the number of parameters (Ai A)» +++ Ap) chosen to be 
representative of the processes involved in a particular physical phe- 
nomenon. Let m represent the number of fundamental dimensions (length, 
time, mass etc.) present in all of the parameters. Then follow these 
steps. 

l. Determine the number of " Tf dimensionless ratios" required to 

incorporate all the chosen parameters (given by n-m). 


2. Combine m+1 parameters into an equation of the following form: 


WTA 
= A A A. 00 (8) 


4 ۱ 


where (A, Age eee An) must together include all of the fundamental 
dimensions. 

3. Substitute the fundamental dimensions into equation (8) and 
equate like powers on each of the fundamental dimensions to zero, 


da Hh Te 
A, Aa rinne mbl = pa 1" t 


4. Solve the resulting m equations with m unknowns to obtain X, , 
y; and 21 ۰ 


5. Substitute the values of Xpypz into the equation for Tf, to 


1 

obtain the first dimensionless ratio. 

6. Repeat steps 2 through 5 after changing A z to another parameter 
m 


not previously used for any A and solve for T. . Repeat as neces- 


sary until T,, Tạ, ۰۰۰۰1۲ و‎ ratios are obtained. 


7. The ratios represented by ۲ T, e. Ton are ail dimensionless 


L? 
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and can be equated one to the other to obtain a functional relation- 
ship that is dimensionally correct. The ratios are then made numer- 
ically equivalent by use of a dimensionless coefficient P, 


An example illustrates the application of the foregoing step by step 


procedure. 
Parameters Dimensions 

Q (total heat present) LT 

W (wind speed) L/t 

fL (coriolis effect) 1/t 

fa (coefficient of thermal expansion) 1/T 

H (mixed-layer depth or MLD) L 

n = 5 m = 3 

STEP 1 n-m = 2 " Mratios” will be required 
STEP 2 T = aw n م‎ and ala aoe li 


To find the powers (x,y,z) which make the above products dimension- 


less: 
STEP 3 
ror m: (L E 

STEP 4 (L) hta o 

(T) *,-1 =0 

(t) zu -% =° 

LEE 

wo ci 

- 1 


2 
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STEP 5 ۱ u, 
r= Qw WA 


OT 
AVE 
im = -z 
+ (9; à 3 
STEP 6 For TE: (im) ۷۶۰۳ 
(L) 4¿ + 2+1 =O 
(T) mh, = 0 
Ct) - یرو وم‎ = 
0 
ka 
Ya = 1 
Ta = q 2w ہر نہ'‎ 
s = = 


STEP 7 Since both um and و1‎ are dimensionless, either can be 


inverted in order to give a relation that agrees with what is observed 


physically. In this case it is known that the MLD is an increasing func- 


tion of the wind, and im and Ta may be combined into: 


1 
TÍ, = 9 
where P is a dimensionless coefficient whose form must be determined 


experimentally. Substituting for m and 1% and solving for the mixed 


layer depth gives: 


MLD = P(Q Bae) 
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In his paper on TI theory, Buckingham [2] gave several methods to 
determine the form of P. Of these, the same method used by Kitaigorod- 
sky [1] will be given as an example here. Since P may be considered a 
function of all the original parameters, they can be used to determine P 
with experimental data. If P does not turn out to be constant, it can 
be made a function of certain of the parameters and the form of that 
function can be found from the data. Thus P is made a function of Y, 
which, in this example, is made to depend upon all the parameters except 
MLD, the quantity to be predicted. Use of TI theory, as before, yields 


a dimensionless ratio of these parameters: 


N ~ BAD 
W 


Enteríng the eguations for both P and N with the same data and plot- 


ting the corresponding pairs gives P(N) in graphical form. 
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APPENDIX II 


RELATIONSHIP BETWEEN VARIANCE OF THE MLD 
AND THE STRENGTH OF THE THERMOCLINE 


Mazeika [6] used 356 groups of BT data to determine that the MLD 
at OWS Echo had fluctuations with a mean amplitude of 17.85 feet and 
Standard deviation of 10.55 feet about the mean. He attributed these 
Fluctuations to internal waves. 

On the basis of 100 BT observations from OWS Papa for the months of 
September and October of the years 1958 through 1960, the author deter- 
mined a system for estimating the amount of fluctuation of the MLD that 
might be expected at verification time. The temperature difference be- 
tween the sea surface and the temperature of the near-isothermal water 
just below the thermocline was used as a measure of the strength of the 
thermocline. The variance of each observation was computed from a five- 
day running mean. The average variance and the average AT for each two- 
meter increment of mixed layer depth was them computed and plotted as 
pairs in figure 7. Figure Y can be used to estimate the variance that 
can be expected at verification time. 

In order to smooth out the fluctuations of the MLD when an accurate 
verification is desired, it is apparently best to take a series of BT 
observations over a period of at least 12 hours and then use the average 


of all the BTs to verify the forecast. 
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APPENDIX Lil 
METHODS USED FOR DETERMINING VALUES OF PARAMETERS 
WIND (W) 

l. Out of the 72-hour period preceding verification time,but not 

including the 12 hours immediately preceding verification time, 

determine the 24-hour period during which the highest average winds 
prevailed or are forecasted to occur. 

2. From the eight wind reports given for that 24-hour period, 

choose the five highest speeds. The average of these five is the 

"representative maximum wind." 

EXCESS HEAT IN UPPERMOST LAYER (Q) 

When finding values of P and N, the computation of Q was based on 
the average BT for each particular i2-hour period. When using equations 
(5) and (7) to forecast an MLD, the most recent BTs available should be 
used to compute 0. In either case the following steps are applicable in 
determining the excess heat in the uppermost layer. 

l. Assume a transitional warm layer exists. For this situation 

compute Qy as follows: 

a) Determine ۸۵۲ و‎ depicted in figure 8. The dashed line 
ís a vertical 0 from the point of maximum curvature of the 
BT trace (roughly the bottom of the thermocline) to the surface. 


b) Qr is then given by 
Q, =F SAREALX 197 (Kg cel/cm”) 


For OWS Papa in summer Poy 5 0.975 for salinity 32.5 °, 00, 


assumed to be constant (Ccp in -cat | and AREA, in MOC). 
OC cm 
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Zi Isa > 3.30% 1074 using Q,, then a seasonal-type thermocline is 
assumed to exist. The excess heat for a seasonal thermocline (04) 
is found by the following objective method, 
a) Determine "AREA." » depicted in figure 9(d). The dotted 
line, figure 9(a), is a vertical drawn from the intersection of 
the BT trace and 200 meters. HH is determined by equalizing 
areas A and B as in figure 9(b). Then 
AREA, 2 CE 


b) EN 
Qe = 7 Sp AREA, X 10 CKg cal /« m? ) 
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Figure 9 


Representation of Area, Used in Determining Qs 
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TABLE 6 


VALUES FOR THE COEFFICIENT OF THERMAL 
EXPANSION BASED ON A SALINITY OF 32.59/oo 


Temperature (°C) 6 < 10° (©) 

> 14.0 205 
14.0 - 13.6 200 
1375 : 111 195 
13.0 - 12.6 190 
12.5 - 1271 186 
12.0 - 11.6 182 
11.5 <, LE 178 
11.0 - 10.6 | 174 
10.5 - 10.1 168 
10.0 - 9.6 163 
9,5 - 9.1 158 
9.0 - 8.6 152 
8.5 - 8.1 147 
8.0 - 7.6 141 

102272) 135.5 

7.0 - 6.0 124.5 
< 6.0 113 
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APPENDIX IV 
FORECASTING METHOD PROPOSED BY KLTAIGORODSKY 

Kitaigorodsky ul, using similarity theory as described by Monin and 
Obukhov [14], developed an equation for forecasting the MLD. In develop- 
ing the equation Kitaígorodsky assumed that thermal convection was neg- 
ligible and that the vertical gradient of salinity was zero. This Limits 
the use of the equation to the warm season when a stable layer exists. 
The equation does incorporate parameters that represent all the processes 
and variables that were considered im section 3 to cause or affect the 
MLD at OWS Papa. 


Equation developed by Kitaigorodsky using similarity theory: 


q 2 
c De _ 
MLD =P SE (8) 
where: Ta = a parameter proportional to the tangentia! 


stress of the wind om the sea surface 


Qp = rate of heat flow from the atmosphere (LT/t); 

g = acceleration of gravity (L/t?); 

coefficient of thermal expansion of sea water (1/T);‏ = ور 
P = a dimensionless coefficient that is a function of‏ 


N; both N and P are determined empirically; 


M > Bl (N is dimensionless); 
65 23 


fL z coriolis parameter (1/t). 
The value of P in equation (8) is given by a plot of 14 correspond- 


ing values of P and N, obtained by Kitaigorodsky from data from the 
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NORPAC ATLAS 023 The large areal extent which the data represents 
makes the results widely applicable throughout the ocean during the warm 
season. 

The present author, using data from OWS Papa [16] has determined 11 
more pairs of values of P and N. Each of these 11 pairs represents from 
four to six BTs averaged together in an attempt to smooth out internal 
wave effects. The paired values of P and N given by both Kitaigorodsky 
and the present author are plotted on figure 10. 

In calculating P and N, it appears on the basis of the reference [17] 
that Kitaigorodsky used an average (climatological) heat flow to find Qp- 
The present author also used an average heat flow, the average heat flow 
since the time that the upper layer of the ocean was last isothermal. 
This average heat flow was arrived at in three steps. 

1. Calculate the heat added in the surface layer since warming 

began. The method described for Q, in Apendix III was used for 

this calculation. 

2. Determine the length of time elapsed since warming began. 

3. Divide value determined in 1 above by the value determined in 

2 above. 

Despite the different locations of the NORPAC and OWS Papa data and 
despite the differences in Qp, both sets of the points P versus N appear 
to fit the same curve (in figure 10). Still both methods for obtaíning 
Qe seem weak. In section 4, Op is replaced by another parameter to 
represent stability in order to eliminate that weakness. 

Another problem in applying the method proposed by Kitaigorodsky 


dP 
lies in the steepness of the slope, qy, for the small values of N. 
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‘Figure 10 | 
Vaiues of P versus N Computed Using Equations Developed by Kitaigorodsky 
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To avoid these problems, the author derived by similarity methods 
an equation with more readily-available parameters and a more stable 


P(N). 
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